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ABSTRACT 


The encounter of the spacecraft Voyager 2 with Neptune and its large satellite Triton 
in August 1989 will provide a crucial test of ideas regarding the origin and chemical 
composition of the outer solar system. In this pre-encounter publication we quantify the 
possibility that Triton is a captured moon which, like Pluto and Charon, originally 
condensed as a major planetesimal within the gas ring that was shed by the contracting 
protosolar cloud at Neptune's orbit. Ideas of supersonic convective turbulence are used to 
compute the gas pressure, temperature and rate of catalytic synthesis of CH 4 , CO 2 and C(s) 
within the protosolar cloud, assuming that all C is initially present as CO. The calculations 
lead to a unique composition for Triton, Pluto, and Charon: each body consists of, by mass, 
18^% solid CO 2 ice, 4% graphite, |% CH 4 ice, 29% methanated water ice and 48% of 
anhydrous rock. This mix has a density consistent with that of the Pluto-Charon system 
and yields a predicted mean density for Triton of 2.20 + 0.05 g/cm^, for satellite radius equal 


to 1,750 km. 
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This report was written at the Department of Mathematics, Monash University, Clayton, 
Victoria 3168, Australia. It was originally submitted for publication to Nature (London) on 
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The encounter of the spacecraft Voyager 2 with Neptune and its large satellite 
Triton in August 1989 will provide a crucial test of ideas regarding the origin and chemical 
composition of the outer solar system. In this report we quantify the possibility 1 ’ 2 that 
Triton is a captured moon which, like Pluto and Charon, originally condensed as a major 
planetesimal within the gas ring that was shed by the contracting protosolar cloud at 
Neptune's orbit 3 ’ 5 . Ideas of supersonic convective turbulence are used to compute the gas 
pressure, temperature and rate of catalytic synthesis of CH^, CO^ and C(s) within the 
protosolar cloud, assuming that all C is initially present as CO. The calculations lead to a 
unique composition for Triton, Pluto, and Charon: each body consists of, by mass, 18-% 
solid CO ice, 4% graphite, -% CH ice, 29% methanated water ice and 48% of anhydrous 

2 2 4 

rock. This mix has a density consistent with that of the Pluto— Charon system and yields a 
predicted mean density for Triton of 2.20 ± 0.05 g/cm 3 . 

The greatest clue to Triton's origin lies in its orbit about Neptune: the orbital 

plane is inclined ~ 159° relative to Neptune's equator, meaning that Triton actually moves 
in a retrograde sense about the planet. Such an irregularity argues against Triton being a 
naturally formed moon of Neptune (N) like the regular satellite systems of Jupiter (J), 
Saturn (S) and Uranus (U) 6 ’ 8 : all of these latter systems orbit the parent body in a 

prograde sense. Instead, it seems dynamically more likely* that Triton was captured from 
a solar orbit after having first condensed as a major planetesimal within the protosolar 
cloud. We shall not pursue the capture process in more detail here, beyond pointing out 
that gas drag within Neptune's primitive envelope of mass ~1M 0 (Earth mass) no doubt 
played a major role in securing this event 9 . Triton's current orbital radius ~ 14.2 R N 
= Neptune's equatorial radius) lies well inside the initial radius ~ 25 R^ expected for 
the proto— planetary envelope, judging from the radial extent of the other regular satellite 
systems. 

Consider now in detail the processes which led to Triton's condensation within the 
protosolar cloud. Many of the basic physical and chemical features of the planetary system 
and regular satellite systems can be accounted for if there had existed a large, radial, 
turbulent kinetic stress within each of the proto— solar/planetary clouds which 


2 


formed these systems 3 ’ 5 ’ 7 ’ 8 . In the simplest representation, </> t v t 2 > = 0pG M(r)/r where p 
is the local gas density, M(r) the total mass interior to radius r, G the gravitation constant, 
and /?= 0.107 is the turbulence parameter 10 ’ 11 . The total radial stress at any point in the 
cloud is then 


p tot = Pgas + <p t v t 2> ’ 

where p gag = pSCljp is the usual gas pressure, T the temperature, p the molecular weight, 
and 9L is the gas constant. Typically, <p^ 2 > ~ 30 p ^ near the outskirts of the cloud, 
meaning that the convective motions are strongly supersonic there. The inclusion of 
turbulent stress has two very important consequences. First, on the physical side, it causes 
a rotating and gravitationally contracting cloud to rid its excess spin angular momentum 
by shedding a discrete family of orbiting gas rings. The orbital radii R n (n = 1,2,3,...) of 
these gas rings, of nearly equal mass m, are related to the mass M c) and 
moment— of— inertia coefficient f (~ 0.01) of the parent cloud by the equations 

V R n+l = t 1 + m /M d f] 2 . (1) 

If the contraction takes place homologously, meaning that both m/M c | and f are constant, 
then the radii R n form a geometric sequence, similar to that of the observed 
planet /satellite orbital radii. 

Second, because the gas pressure in the outer layers of the turbulent cloud 
(P — 0.107) is very much less than that in conventional non— turbulent models, carbon 
chemistry is shifted into a region where the formation of C0 2 and solid carbon C(s) 
becomes possible. At low temperatures, the carbon activity a c is controlled by the 
reaction CH ^ 2H + C(s). We have 

4 2 w 

a C = [ p CH / p H 2 1 K ,( T ) > ( 2 ) 

4 2 

where the pds denote partial pressures and K^T) is the equilibrium constant 12 . For a fixed 

ratio of CH^ to H^, a^, varies as 1 /Pjj • In & gas of solar composition 13 , graphite can form 

2 

( a C - 1) onl y if T < 473 K and p^ <4x 10" 8 bar. At 470 K, a^ for the conventional 
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solar adiabat 12 is 5 * 10' 4 , since p H ~ 2 * 10' 4 bar is so large. That is, graphite formation 

2 

is impossible here. 

Within a turbulent cloud of solar composition whose equatorial radius matches that 

of Neptune’s orbit, a c (470 K) = 0.82, when p H = 7.0 * 10~ 8 bar. Graphite formation is 

2 

thus almost possible. For such a cloud, however, having heavy element abundance Z = Z 0 , 
the total mass of rock and ice condensate within the gas rings shed at the orbits of J, S, U, 
and N are 9.5, 7.1, 5.2, and 3.8 M 0 , respectively. If the condensate in each ring were to 
aggregate into a single planetary core, the masses of the cores would fall short of the values 
measured by Pioneer 11 (ref. 14), or deduced from planetary interior models 15 , by factors of 
2.1, 2.6, 2.5, and 4.2, respectively. The same problem occurs in the Galilean system of 
satellites, where eqn. (1) yields a mass for Io of 3.4 * 10 25 g. This falls short of the 
observed mass 16 by a factor of 2.6. To remedy this deficiency, which is probably due to the 
simplified treatment of supersonic turbulent stress, it was decided to increase Z by a factor 
of 2.5. Under these circumstances a c > 1 for all cloud equatorial radii R e > 3,300 R @ , 
provided thermodynamic equilibrium is achieved. 

Consider now the rate of production of CH^, C(s) and CO^ in a cloud which is 
initially composed of CO. The controlling reaction is CO + 3H^ ^ CH^ + H^O. For 
simplicity, we assume a linear equation similar to the Temkin— Pyzhev rate equation for 
iron— catalysed NH^ synthesis 17 ’ 18 : 

dp CH /dt = k [p C o - P CH P H o/ p H 3 K 2 ^ ‘ 

4 L 4 2 2 


Here k , K are the rate and equilibrium constants. Collision theory 19 yields 

2 2 


k 

2 


3/?X Fe,f 

/, Fe a Fe 


#T 

^CO 




E a) 


( 3 ) 


where Xp g ^ = 0.487 Xp g is the mass fraction of free metallic iron (allowing for FeS 
formation), pp e = 7.875 g/cm 3 , a Fg is the iron grain radius, and E a = 21.3 kcal/mole is 
the activation energy 20 . Although earlier calculations 21 ’ 22 assume that a Fe = 10‘ 2 cm — a 
mean derived from studies of chondritic meteorites 23 — we feel that a much smaller value 
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is likely to apply to the hot convective interior of the protosolar cloud. Unless 
a Fe ~ ^ 5 cm ’ rate ^ P r °duction is negligible and, if ap e >2* 10' 6 cm, we 

cannot account for any free CH ice on Pluto's surface 24 because H O enclathrates all the 

4 2 

available CH^. We therefore adopt ap g = 10" 6 cm. 

Eqn. (3) can be integrated once the radial velocity v g = dR g /dt and initial radius 
R^ of the cloud are specified. Now for R 0 < R* ~ 426 R @ , the cloud's gravitational energy 
&grav - 3-3 GM C 2 /R e exceeds the total energy # e q m needed to establish a complete 
hydrostatic and thermodynamic equilibrium. Radial contraction is governed by the rate at 
which the excess energy A# = _ can be radiated from the surface. Very 

roughly, d(A#)/dt ~ - (GM c 2 /R 2 )Vg = where 2.82jra s R 2 Tgf£ is the luminosity 11 , 
and ~ 361 (R*/R e ) 0 ' 9 K is the effective surface temperature. We obtain 

v e ~ 0.0103 [R e /R*] 0 ' 4 Rgj/yr . (4) 


We assume that eqn. (4) is also valid for R g > R*, noting that the formation of a small 
compact core at the centre can release sufficient energy to stabilize the rest of the cloud 3 . 

Fig. 1 shows the run of methane number fraction f^ for a cloud whose initial 


4 

radius is = 7946 R 0 . This radius was chosen so that the composition of the condensate 
in the gas ring shed at Neptune's orbit has a compressed density consistent with the 
observed mean density of the Pluto— Charon system (see below). Choosing R^ > 7946 R 0 
yields a condensate which has too much CH ice and too low a density. Now f^ is taken 

4 

to equal the maximum CH^ fraction within the cloud of given size Rg. At the cloud 
surface, where T is least, the rate of production of CH is negligible. In the deep interior, 

4 

where the rate is large, the thermodynamically favoured state for C is CO. That is, there 

is a unique depth and temperature T where f^u is a maximum. The values of T 

max kj ri max 


4 

for cloud sizes R g which match the orbits of J, S, U, and N are 458, 465, 470 and 484 K, 

respectively. We assume that convective mixing homogenizes the outer layers of the cloud, 

yielding a uniform composition having mean CH fraction T^t, = f^rr . Almost 

4 vj'Ii L/ii .max 

4 4 

certainly, some CH will be destroyed through downward mixing into regions where CO is 

4 
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dominant. That is, is probably less than m ax' ^ more comprehensive study of 

4 4 ,m 

the effect of mixing on the CH level is outside the scope of this report. 

4 

The rate of production of solid carbon C(s) is assumed to be proportional to that of 
CH 4 - Thls leads t0 f C(s) = f CH t f C(s)/ f CH leqm’ where f C(s),eqm’ f CH ,eqm denote the 

4 4 4 

corresponding equilibrium fractions of C(s) and CH at temperature T . Lastly, the 

4 max 

ratio f^Q / f^Q is assumed to equal its equilibrium ratio. This is valid owing to the fairly 
2 

rapid rate of the gas phase reaction CO + HO^ CO + H (ref. 25). 

2 2 2 

At the orbit of Neptune f CH “ °- 153, f C(s) = 0111 and f CO = ° 143 ' These 

4 2 

fractions enable us to calculate the abundances of anhydrous rock 
(MgSiO ,Fe 0 ,FeS,Al 0 ,...), H O, etc., in the Neptunian gas ring, noting that O, which 
is normally tied up as H 0 in a CH —rich atmosphere, or between H 0 and CO in a 

2 4 2 

CO— rich atmosphere 21 ’ 26 , is now to be first apportioned to CO and CO^. Assuming the 
Anders and Ebihara heavy element abundances 13 , multiplied by the factor 2.5, the 
percentage masses of rock, ^0, CH^ C0 2 and C(s) are 1.321, 0.677, 0.201, 0.507 and 
0.109%, respectively. All of these materials, except CH , fully condense from the gas 

4 

phase, since the gas ring temperature T^j = 27.4 K is so low. Of the CH^, 52.2% is 
encaged by H O ice, to form CH »5.75 H O, whilst only a modest 6.3% can condense out as 

2 4 2 

solid CH^ ice since the condensation temperature only just exceeds Tj^ (by 0.1 K!). 

The mass fractions of the condensate at Neptune's orbit are thus: rock 0.4835, H^O 

ice 0.2478, clathrated CH^ 0.0384, CH^ ice 0.0046, C0 2 ice 0.1857 and graphite 0.0400. 

This is the predicted bulk composition of Pluto, Charon, and Triton, assuming that all 3 

objects condensed from the Neptunian gas ring. No NH ice is present since the rate of 

3 

NH synthesis in the protosolar cloud is negligible. That is, all nitrogen exists as N . Since 
no room remains for the incorporation of CO, N 2 and Ar in the H 2 0 ice, after the 
energetically favoured 27 enclathration of CH has taken place, the satellites should contain 

4 

very little trace of these species. That is, the atmosphere of Triton should consist mostly 
of CH . 

4 
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Structural and evolutionary models may now be constructed to estimate the 

satellite mean densities p ga ^. These calculations are described only briefly here. First, a 

set of cold chemically homogeneous models for Charon and Pluto were constructed having 

radii of 596 and 1142 km, equal to the observed values 28 . The assumed uniform 

temperature is T g = 45 K (ref. 29). The values Pq^ = T9984, Ppj = 2.0744 g/cm so 

obtained have a mass-weighted average of 2.0649 g/cm 3 , which coincides with the observed 

system mean p g ^ g = 2.065 ± 0.047 g/cm 3 . This justifies the choice R.^ 7946 Rg, 

mentioned earlier. Next, the cold satellite models were allowed to warm up under the 

influence of radiogenic heating and thermally evolved for a period of 4,600 Myr. This 

yielded a present-day temperature profile whose central temperature T , along with the 

central pressure P , is shown in Table 1. The corresponding values p gat are also shown. 

The densities are only slightly less than those of the cold, isothermal models. By choosing 

a slightly smaller initial cloud radius R q 2 7890 R @ , one could achieve a slightly less icy, 

and hence denser, model. An estimate of the correct mean densities 

p =p + 0.0066 g/cm, whose system average coincides with the observed value, is 

shown in the last column. It should be mentioned that the model of Pluto reported here is 

cooler than that of McKinnon and Mueller 2 . This is partly due to a lower rock fraction 

(0.48 vs 0.73), which reduces chondritic heating, and partly because of a higher thermal 

conductivity, brought about mostly by the presence of graphite 30 . 

The final 3 rows of Table 1 list the computed properties of 3 different-sized models 

of Triton. The radii span the observed range r^ r = 1,750 ± 250 km 31 . The surface 

temperature is 57 K (ref. 32). Triton is so large that, in the absence of internal convection, 

T rises to 670 K after 2,800 Myr and equals 581 K at present age. That is, wide— spread 
0 

melting of ice and internal differentiation of rock and ices would ensue, resulting in a 

substantial drop in mean density (from 2.2 to 1.8 g/cm 3 ). We suggest here that the 

presence of C0 2 ice prevents this outcome. CC> 2 ice is a weakly bonded crystal whose 

isothermal compressibility 33 Xqq at 150 K is some 2^ times that of f^O ice, for the 

2 

pressures applicable to Triton's interior. The ratio X^q lx jj q increases with T. We 

2 2 

therefore suggest, following Lupo and Lewis' argument for weak CH 4 ice 34 , that solid-state 
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convection quickly wipes out any steep conductive temperature gradients inside Triton 

35 

and, thereafter, controls the thermal evolution. To model this phenomenon we assume 
that creep sets in as soon as T rises to 0.6 of the CC> 2 melting temperature T m (ref. 36). 

On this basis, Triton has a central convective core covering ~ 65% of its mass. The 
temperature T , T at the centre and core edge are given in the Table, along with /> gat 
and the corrected value p corr . From these results we make the prediction that, if Voyager 
2 finds Triton's true radius to be 1,750 + Ar km, the observed mean density p Tr will satisfy 
the equation 

p Tr = 2.20 + 3.00 x 10' 4 (Ar) + 2.0 x 10‘ 7 (Ar) 2 ± 0.05 g/cm 3 . 

The author thanks J.D. Anderson, J.R. Anderson, A.F. Moody and G.K. White for 
valuable discussions, and Telecom (Australia) for financial support. 
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Table 1 Structural properties of COg-rich model satellites 


r sat 

P 

0 

T 

0 

T core 

m core/ m sat 

^sat 

p 

r corr 

(km) 

(kbar) 

(K) 

(K) 


(g/cm 3 ) 

(g/cm 3 ) 

596 

2.05 

61 

— 

0 . 

1.993 

1.999 

1142 

8.05 

132 

— 

0 . 

2.068 

2.074 

1500 

15.5 

239 

181 

0.494 

2.133 

2.140 

1750 

23.0 

271 

180 

0.649 

2.196 

2.203 

2000 

33.6 

308 

181 

0.739 

2.284 

2.290 
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Figure 1. Relative number fractions f, in terms of total carbon, of CH 4 , CO, COj and C(s), in 
the outer layers of the turbulent protosolar cloud, plotted against cloud radius Rg, 
measured in units of Rq. C(s) means graphite. Rq is the initial cloud radius, (J 

the turbulence parameter 5 ' 10 , Z the total heavy element fraction, and ap e is the 
mean radius of the iron grains used in the calculation of the rate of catalytic 
production of CH 4 from the initial reservoir of CO. 








